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Abstract 

We report a search for muon neutrino disappearance in 
the Am 2 region of 0.5 — 40 eV 2 using data from both Sci- 
BooNE and MiniBooNE experiments. SciBooNE data 
provides a constraint on the neutrino flux, so that the 
sensitivity to disappearance with both detectors is 
better than with just MiniBooNE alone. The prelim- 
inary sensitivity for a joint disappearance search is 
presented. 

1. Introduction 

Neutrino oscillations have been observed and con- 
firmed at mass splitting (Am 2 ) of - 10~ 5 eV 2 and 
~ 10~ 3 eV 2 , called the "solar" and "atmospheric" re- 
gions, respectively. The observed mixing is consistent 
with three generations of neutrinos. 

However, the LSND experiment observed an excess of 
v e in a v ^ beam, indicating a possible oscillation in the 
Am 2 ~ 1 eV 2 region [1] . To explain LSND with oscilla- 
tions requires more than three generations of neutrinos 
or other exotic physics beyond the Standard Model. 

To test the oscillation at Am 2 ~ 1 eV 2 , the Mini- 
BooNE experiment recently made searches for both v e 
appearance [21 [3] and disappearance [4] in this pa- 
rameter region. The experiment observed no signifi- 
cant v e appearance signal and ruled out as being due 
to 2-neutrino oscillations. However, the sensitivity of 
MiniBooNE-only disappearance search was limited 
by the large flux and neutrino interaction cross-section 
uncertainties. 

Here, we discuss an improved search for disappear- 
ance using data from both the SciBooNE [5] and the 
MiniBooNE experiments, where SciBooNE detector is 
used to constrain flux and cross-section uncertainties. 

2. Experimental Setup 

2.1. Fermilab Booster Neutrino Beam 




100 m 440 m 

Fig. 1. The setup of SciBooNE and MiniBooNE experiments. 

The experiments use the Booster Neutrino Beam 
(BNB) at Fermilab [6]. The primary proton beam, ex- 
tracted with a kinetic energy of 8 GcV, strikes a 71 cm 
long, 1 cm diameter beryllium target. The mesons, pri- 
marily 7r + , generated by the p— Be interactions are fo- 
cused with a magnetic horn and decay in the following 
50 m decay volume, producing an intense neutrino beam 
with the peak energy of ~0.7 GeV. When the horn po- 
larity is reversed, tt~ are focused and hence a predomi- 
nantly antincutrino beam is created. 



2.2. SciBooNE Detector 

The SciBooNE detector [5] is located 100 m down- 
stream from the beryllium target. 

The detector complex consists of three sub-detectors: 
a fully active fine grained scintillator tracking detec- 
tor (SciBar), an electromagnetic calorimeter (EC) and 
a muon range detector (MRD). 

The SciBar detector consists of 14,336 extruded plastic 
scintillator strips (CH), each with dimension of 1.3 x 
2.5 x 300 cm 3 . The scintillators are arranged vertically 
and horizontally to construct a 3 x 3 x 1.7 m 3 detector. 
The detector itself is the neutrino target and its fiducial 
volume is 10.6 tons. 

The EC is installed downstream of the SciBar, and is 
made of scintillating fibers embedded in lead foil. 

The MRD is located downstream of the EC in order to 
measure the momentum of muons up to 1.2 GeV/c using 
the muon range. It consists of 12 layers of 2"-thick iron 
plates sandwiched between layers of 6 mm-thick plastic 
scintillator planes. 

The SciBooNE experiment ran from June 2007 until 
August 2008, collecting a total of 2.52 x 10 20 Protons 
on Target (POT) for physics analysis; 0.99 x 10 20 POT 
in neutrino mode and 1.53 x 10 20 POT in antincutrino 
mode. 



2.3. MiniBooNE Detector 

The MiniBooNE detector [7] is located 440 m down- 
stream from the SciBooNE detector. The detector is a 12 
m diameter spherical tank filled with 800 tons of min- 
eral oil (CH2). The MiniBooNE experiment has been 
taking beam data since 2002, including the SciBooNE 
and MiniBooNE joint-run period. The collected number 
of POT after data quality cut in the neutrino mode is 
5.579 x 10 20 in addition to the data from the joint-run 
period. 

3. Disappearance Analysis 

3.1. Analysis Overview 

In this paper, we report only the neutrino data [y^ — > 
u x ) disappearance analysis. We search for muon neu- 
trino disappearance by comparing neutrino fluxes at Sci- 
BooNE and MiniBooNE detectors. 

The analysis is performed in the following three steps: 
(1) Neutrino flux measurement at SciBooNE, (2) Flux 
extrapolation to MiniBooNE, and (3) Oscillation fit. 

At each step, systematic errors are estimated and 
propagated to the final prediction. The majority of the 
flux and cross-section uncertainties cancels since the neu- 
trino interaction target in both detectors is effectively 
carbon, and the two detectors are on the same beam 
line. 

We describe these steps in detail in the following sec- 
tions. 
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3.2. Neutrino Flux Measurement at SciBooNE 

Charged Current Event Selection 

For the spectrum analysis at SciBooNE, we use inclu- 
sive charged current (CC) interactions, whose signa- 
ture is long muon tracks. First, we identify muons by se- 
lecting the longest track with energy deposit consistent 
with a minimum-ionizing particle. Second, we require 
the vertex of the track to be within the SciBar fiducial 
volume. The events are further divided into two sub- 
samples based on the location of the muon track end 
points: a "SciBar-stopped" sample containing muons 
that have stopped inside the SciBar detector and a 
"MRD-stopped" sample with muons that have stopped 
in the MRD. These two samples each contain approxi- 
mately 14k and 20k events with mean energies of 0.8 and 
1.1 GeV, respectively. 



Spectrum Fitting 

The neutrino spectrum at SciBooNE is extracted by 
fitting muon momentum (P^) and muon angle (9^) dis- 
tributions from each sample. 

We prepare MC templates for P^ and <9 M distributions 
for several true neutrino energy (P„) regions. The E v 
regions are divided by 250 McV up to 1.25 GeV, and a 
single region is assigned for E v > 1.25 GeV. Then, the 
scale factors for each E v region are determined to mini- 
mize the x 2 between data and MC. Figure [2~H shows the 
fit result. The systematic errors from SciBooNE detector 
response and neutrino cross-section models are estimated 
and shown in the plot. 

Figure 151 is the P M and 9^ distributions of SciBooNE's 
MRD-stopped sample, after applying scale factors ob- 
tained by the spectrum fitting. We confirm the MC dis- 
tributions agrees well to data after fitting. 
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Fig. 2. Scale factors obtained by SciBooNE spectrum fitting. 
The error bars show the sum of SciBooNE statistical and 
systematic uncertainties. 



3.3. Flux Extrapolation to MiniBooNE 

MiniBooNE Event Selection 

We select events in MiniBooNE by requiring single 
muon and its decay electron. Neutrino energy is recon- 
structed from muon kinematics by assuming CC Quasi 



Elastic (CCQE) interaction {v^n —> fi p): 
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Fig. 3. Distribution of reconstructed muon momentum (top) 
and muon angle (bottom) for the MRD-stopped sample. The 
dots show the data, and histograms show the MC prediction 
with the contributions from neutrino interaction modes. The 
MC distributions are tuned by the E„ scale factors obtained 
by the spectrum fit. 



where AM = M 2 — M 2 ; M indicate the muon, proton, 

or neutron mass with appropriate subscripts; Eb is the 
nucleon binding energy; E^ is the reconstructed muon 
energy. 

MiniBooNE E^ ec prediction 

To predict the E^ ec distribution at MiniBooNE, we 
extrapolate the measured SciBooNE flux to MiniBooNE 
in two steps. 

First, we apply MiniBooNE/SciBooNE flux ratio to 
make a prediction of the true neutrino energy distribu- 
tion at MiniBooNE. Then, we smear the true neutrino 
energy prediction to the reconstructed neutrino energy. 

Systematic uncertainties for the flux ratio is estimated 
by varying the cross-section and flux models. Addition- 
ally, the uncertainties of the smearing function, which 
convert true E u to E^ ec , is estimated by varying the 
cross-section models. 

Finally, we add MiniBooNE detector response error to 
the E^} ec prediction. 

The predicted MiniBooNE reconstructed neutrino en- 
ergy distribution and its systematic uncertainties are 
shown in the Figure |4~1 



3.4- Oscillation Fit and Sensitivity 
Fit Method 

We test the oscillation hypothesis assuming the mixing 



between 2 neutrino flavors; 



and 



The 



disappearance probability is given as 

PK -» Vx) =sin 2 26>sin 2 (1.27ATO 2 J L/P), 

where 9 is the mixing angle, Am 2 [eV 2 ] is the mass split- 
ting between 2 flavors, L[km] is the distance traveled and 
E[GcV] is the neutrino energy. 
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Fig. 4. Predicted MiniBooNE reconstructed neutrino energy 
distribution. MiniBooNE detector error, flux and cross-sec- 
tion uncertainty, and the total systematic uncertainty are sep- 
arately shown. 



We fit the MiniBooNE E* ec distribution to find the 
best fit parameter minimizing the \ 2 value: 



16 bills 



= J2 (N? ata — Nf)M^ 1 (N^ 



da t o 



where i,j denote E^ ec bins, Nf°j ta and Nf- denote ob- 
served and predicted number of events at each bin, re- 
spectively, and My represents statistical and systematic 
uncertainties for the final E^ ec prediction. 

Then we define the allowed region by A\ 2 = 
X 2 {true) — x 2 (best) values, where \ 2 (true) is the % 2 at 
the oscillation prediction being tested, and % 2 (6es£) is 

the smallest \ 2 value across the (Am 2 , sin 2 29) plane. 

To obtain the confidence level at each oscillation pa- 
rameter point (Am 2 , sin 2 28), we use Feldman- Cousins' 
method jS]. In this method, 1000 "fake-data" predic- 
tions are formed, using random draws of the statistical 
and systematic uncertainties and some underlying oscil- 
lation hypothesis. Then, each fake-data is fit to obtain 
the relation between the A% 2 values and the correspond- 
ing probabilities. This process is repeated for each pair 
of (Am 2 , sin 2 26) true oscillation parameter being tested. 



Expected Limit 

The sensitivity is defined as the average of limits ob- 
tained from fake experiments with null underlying oscil- 
lation. 

Figure I5~1 shows the 90% CL. sensitivity for the dis- 
appearance. The expected ±1<7 band is also shown in 
the plot. The expected sensitivity directly supersedes 
the MiniBooNE only disappearance result, as sub- 
stantial flux and cross section uncertainties have been 
reduced. 

4. Summary and Prospects 

We present SciBooNE-MiniBooNE joint analysis of a 
search for disappearance in a accelerator neutrino 
beam. The analysis is sensitive to the oscillation at the 
Am 2 region of 0.5 — 40 cV 2 . The sensitivity to disap- 
pearance has been improved relative to the MiniBooNE 
shape-only analysis, with results to be released soon. In 
addition, a joint anti-neutrino oscillation analysis will be 
performed using the anti- neutrino data set. 
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Fig. 5. The expected sensitivity for f M disappearance. The 
dotted curve shows the 90% CL limits from CDHS [9] and 
CCFR [10] experiments. The thin solid curve is the Mini- 
BooNE-only 90% CL sensitivity. The thick solid curve and 
the filled region are the 90% CL sensitivity and ±lcr band 
from SciBooNE-MiniBooNE joint analysis, respectively. 
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